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Gravitational instability of expanding shells 
Solution with nonlinear terms 
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Abstract. A model of the thin shell expanding into a uniform ambient medium is developed. Density perturbations 
are described using equations with linear and quadratic terms, and the linear and the nonlinear solutions are 
compared. We follow the time evolution of the fragmentation process and separate the well defined fragments. 
Their mass spectrum is compared to observations and we also estimate their formation time. 
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1. Introduction 

Hq, emission, which is frequently observed along the pe- 
riphery of giant expanding HI shells, is related to star for- 
mation, which has been triggered due to the compression 
of the ISM. As an example of correlated star formation, 
the ring of young stars observed near to the Sun consist- 
ing of Orion, Perseus, Sco-Cen and other OB associations 
can be used (Lindblad et al., 1997). The HI hole No. 35 
in the dwarf galaxy IC 2574 (Waher & Brinks, 1999) with 
the Hq emission found along its rims serves as another 
example of triggered star formation. 

The analysis of the stability of a pressure-confined slab 



performed by Elmegreen & Elmegreen (1978) has been ex- 
tended to spherical expanding shocks by Vishniac ( 19831 ). 
The time development of the gravitational collapse of lin- 
ear perturbations in decelerating, isothermal shocked lay- 
ers has been examined numerically and analytically by 



Elmegreen (1985, 1994). The purely hydrodynamical non- 
linear instability has been discussed by Vishniac (1994). 
In this paper, we continue with nonlinear analysis of the 
gravitational instability of spherically symmetric shells ex- 
panding into stationary homogeneous medium. 

We modify the approach adopted by Fuchs ( 199(j ) who 
described the fragmentation of uniformly rotating self- 
gravitating disks. If some conditions are fulfilled, the ex- 
panding shell may become gravitationally unstable and 
break to fragments. The inclusion of higher order terms 
helps to determine with better accuracy than the linear 
analysis when, where and how quickly it happens. 

The Rayleigh- Taylor (R-T) instability is not expected 
to develop in the situation explored because a spheri- 



cal shock expanding into the homogeneous interstellar 
medium is always decelerated collecting stationary am- 
bient medium. R-T instabilities may be important in dif- 
ferent situations when the density of the ambient medium 
drops down sufficiently quickly so that the shell can ac- 
celerate mixing hot and cold gas components. This is the 
case of very active SF regions where the shell interacts 
with previously formed fragments. The region 30 Doradus 
in LMC may serve as an example of R-T instability in 
action as described by Redman et al. (1999). 

For highl y sup ersonic flows multiple shocks may de- 
velop (Falle, 



198ll) , but at that time the shells are gravi- 



tationally stable due to squeezing connected to the fast ex- 
pansion. Later, when they decelerate to velocities less than 
50 kms~^, the gravity starts to be important, while radia- 
tive instabiliti es of t he outer shock described by Strickland 
and Blondin ( 1995| ) loose their influence. 

This work may be extended to nonspherical oscil- 
lations using the formalism worked out by Bicak & 
Schmidt ( 1999 ) for cosmological applications. Here we also 
ignore the deviations from spherical symmetry resulting 
from initial asymmetry of the energy input. In a smooth 
medium with only large scale density gradients the shell 
approaches quickly the spherical symmetry as demon- 
strated by Bisnovatyi-Kogan and Blinnikov ( 1982 ). Radial 
perturbations resulting from inhomogcneitics of the am- 
bient medium and variations of the shell surface density 
will be discussed in a subsequent paper. 
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2. The expanding shell in a static, homogeneous 
medium 

The energy input from an OB association or other sources 
creates a blast-wave propagating into the ambient medium 
(Ostriker & McKee 1988; Bisnovatyi-Kogan & SiHch 
1995). Since its radius is during the majority of the evo- 
lution much larger than its thickness, the thin shell ap- 
proximation can be used. The blast- wave is considered as 
an expanding infinitesimally thin layer surrounding the 
hot medium inside. The analytical solution of the expan- 
sion in the thin shell approximation was derived by Sedov 
(1959). In a static, homogeneous medium without the ex- 
ternal or internal gravitational field, the blast-wave is al- 
ways spherically symmetric, it sweeps the ambient matter 
and decelerates. Its evolution is given by the equation of 
motion 

^(AfT/) = 5(Pint-Pcxt), 

the equation of mass conservation 
d 
di 



M = SpoV, 



(1) 



(2) 



and the equation of state 
" 3 T^' 

where M and V are the total mass collected in the shell 
and its expansion velocity, S is the shell surface. Pint and 
Poxt are pressures inside and outside of the shell, po is the 
density of the ambient medium and Eth and Vol are the 
thermal energy and volume inside of the shell. Neglecting 
the external pressure and keeping the Eth — const, and 
assuming Eth = f^'tot, where i?tot is the total input en- 
ergy from the source, the radius of the shell R grows with 
time as 

R\ \i/5 



72.2 
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V 1cm ^ / 

it decelerates as 
V 
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and its unperturbed surface density Sq 
as 



Mq pc- 



0.57 



' PO 

V Icm^^ 
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(5) 

grows with time 



(6) 



The linear analysis of the gravitational instability of 
expanding shells by Elmegreen (1994) gives for the in- 
stantaneous maximum growth rate ojbge of a transversal 
perturbation of a shell 

r OH 1/2 

' R ^ i?2 ^ 



(7) 



where c is the sound speed within the shell. The instability 
occurs for wbge > 0. Inserting Eqs. (^ - ||) into Eq. (^, 
we may derive the time when the instability occurs for 
the first time. WBGE(ib) = for 



VMyr 



28-8 (l^)'/' 



\ IDS' 1 erg J \ Icm-i' / 



(8) 



The ratio of the wavelength A of the fastest perturbation 
to R 



X/R 



2c2 



GPSn 



(9) 



is at t > tb less than -^^^ . The sound speed c in the dense 
and cold shell is always smaller than the speed of sound in 
the ambient medium, which is smaller than the expansion 
speed of the shell V. Therefore, X/R << 1. 

In a subsequent paper, we shall also discuss the insta- 
bility in non-spherical shells: the values of R, V, and Sq 
will be taken from numerical simulations. 

3. Hydrodynamical and Poisson equations on the 
surface of the shell 

We consider the cold and thin shell of radius R surround- 
ing the hot interior and expanding with velocity V into a 
uniform medium of density po . The intrinsic surface den- 
sity of the shell S is composed of unperturbed part Eg 
plus the perturbation Ei (E = Eq -I- Ei). Perturbation 
El results from the fiows on the surface of the shell redis- 
tributing the accumulated mass. We assume that Eq cor- 
responds to i? as Eg = poR/3, which means that all the 
encountered mass is accumulated to the shell. (It comes 
from Eo = ^ttR^ po/^irR^). 

The mass conservation law in a small area on the sur- 
face of the shell is 



dm , . 

— + {V,mv) = AVpo, 



(10) 



where m is mass in the area A = AnaR^, a is a fixed 
small fraction of the sphere, v denotes a two-dimensional 
velocity of surface flows connected to perturbations above 
the stretching of the area due to expansion of the unper- 
turbed shell: r and v — r are two-dimensional vectors 
in the tangential plane of the shell at the central point 
of the considered area. We consider angular coordinates 
= r/R and angular velocity ft — v/R to describe the 
evolution on the surface of the shell (see Fig. 0) . 

With E = m/A we obtain continuity equation in a 
form 

^ + 2E-^-f Ei?(V,f2)-f P(0,V)E-Fpo , (H) 
ot R 

The equation of motion for flows on the shell surface 
has a form 



1 d(mv) 9 

^ ' = -c^VE - EV$, 



A dt 



(12) 
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Fig. 1. The coordinates on the shell surface: © = r/i? for 
the position and the angular velocity O = v/R ioY the 
surface flows. 

where c is the constant isothermal sound speed inside the 
cold shell, $ is the gravitational potential generated by 
the mass distribution given in the tangential plane to the 
shell. The shell is confined by the pressure of the bubble 
interior and the outer shock. We assume that the pres- 
sure distribution through the shell corresponds to constant 
isothermal sound speed c and constant volume density in- 
side of the shell with a slight deviations close to its inner 
edge, which is in contact with the hot interior. Using the 
continuity equation (pT|), we can rewrite the Eq. (n2[) as 



+ vn + wn— + R^ (n, v)n- R^n (v, n) 

at S 



= VE - V$ 



(13) 



The gravitational potential <& is related to the surface 
density by the Poisson equation 

A$ ^ 4TrGY.S{z) , (14) 

where G is the constant of gravitation and S{z) is a delta 
function of the space coordinate z perpendicular to the 
surface of the shell. 



4. Perturbation analysis 

We assume a small perturbation of the shell surface den- 
sity El <C Eq which evolves due to surface flows given with 
velocity v. The perturbed hydrodynamical equations (pT|), 
( p^ ) and the perturbed Poisson equation (^^ have form 

^^^+2E ^ 



Ei? (V, J7) + i? (rj, V) El = 0, 



(15) 



R^ + R^ in, V) n - R^n (V, vt) 

at 

= -^(l-^)VEi-V$i-4yf2-f3yn|i, (16) 



A$i = 47rG'Ei(5(z), (17) 
where the 1/E in Eq. (|l|) was evaluated up to quadratic 



terms in Ei 



rVE 



higher order terms 



) VEi, and VEq = on the shell surface. 

The perturbation of the surface density Ei and the 
angular velocity of the surface flows $1 can be written as 



El 



,i{r,,&) 



(18) 



which may be inserted to the perturbed Eqs. (|15|) and 
(|l6|). T] denotes a dimensionless wave- vector ij = kR. We 
assume no surface macroscopic flow through all the con- 
sidered area which means fio = 0. Further we assume 
5^10 = (mass accumulation due to expansion to the am- 
bient medium is included in Eq). The Fourier transform 
of the Eq. (|l|) is 

E^ +Eo (-ZJ7, n^n) + 2]^E^ + J^r,' ^^-r,') ^r,-r,' 

where we used the identity t]' + [t] — rj') = -q. 
The solution of the Poisson equation dlT 



V$i = -27rGX E^^e'^"^®) (20) 

can be inserted to the Fourier transform of the Eq. (p^. 
We get: 



v' 



3V 



(21) 



5. Linear analysis 

We give the solution of Eqs. (|9|) and (|2l| ) using the Hnear 
terms only. This linear solution will be later compared to 
the results obtained with the nonlinear terms. The com- 
plete linear analysis of the expanding shell was also done 



by Elmegreen ( |1994[ ). 

Linearized Eqs. ( [l^ ) and (21) have a form 



E^ + Eo(-ir7,0_^) + 2-^E^ = 

2 

RCtr, = -—ikY^r, + 27rG^E„ - WVlr,. 
So ' Iryl ' 

Angular velocity rj.^ can be split in two components 
parallel and orthogonal to the wave- vector rj 

n -n ^ + n ^ 



(22) 



(23) 



(24) 
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where r]± is a vector in the tangential plane, which is 
perpendicular to the wave-vector rj and |?7^| = jryj. 
We get the set of equations 



V 



2t:G 



(25) 
(26) 
(27) 

(28) 



We find the eigenvalues and eigenvectors of the linear 
operator C 



which can be formally written as 







( 















V 



y2 ^2^2 27rGI]o?7 



i?2 



R 



V 



(29) 
(30) 



where rj is the magnitude of r). 
The related eigenvectors are 



(1,2) 



(3) 



(31) 



(32) 



oj'^^ indicates decrease, 



and it have meaning of the perturbations growth rate. 

(2) 

The time evolution of the imaginary part of the is 
shown by Fig. p[ 

Since the eigenvalues uori ' depend on rj through the 



,(1,2) 



relation (|29|), the maximum perturbation growth rate 



can be found: 



and occurs at the wavenumber: 

^max — 9 ■ 



(33) 



(34) 



30 40 
t [Myr] 



,(2) 



Fig. 2. The time dependence of the imaginary part of the 
which can cause the instability. The Sedov solution 
was used with following parameters: total energy Etot = 
lO^^erg, density of ambient medium hq — 1 cm~'^, average 
molecular weight /i = 1.3, sound speed in the shell c = 
1 kms"\ 



(2) n K 

77, max is close to ujbge derived by Elmegreen (0994), 



actually Wj^^max — — ici^BGE- If the shell is unstable, i. e. 

(2) 

the imaginary part of the w,, is negative, for the certain 
time, the shell may break to fragments. In Ehlerova et 
al. (1997) is the fragmentation integral for the maximum 
perturbation growth rate defined as: 
t 



(35) 



where th is the time when the instability begins. We gen- 
eralize it for any unstable mode with the wavenumber ry 
using an analogy to the value for wbge as 



The eigenvalues ujrj are time dependent and can be 
used to obtain a criteria for the instability of the shell. For 
a short time for which the change of the linear operator C 
is small (and we can always find the time interval which 
is short enough), the solution of Eqs. ( p8|) is a part of the 
exponential function ~ e''^*. The has always meaning 

f 1 2^ 

of the decrease of perturbations. If Luij have got a real /f,r)(^f,r)) 
part, solution is stable with decreasing oscillations. If not, 

uj^P can be imaginary negative 



-lUJ. 



(2) 



{t')dt' 



(36) 



The fragmentation time tf _,, (the time giving some devel- 
opment level of a fragment) is defined as the time when 
the fragmentation integral is equal to one: 



,(2) 



{t')dt' = 1 . 



(37) 



6. Nonlinear analysis 



The equations with nonlinear terms are so lved by the sim- 
ilar procedure as adopted by Fuchs (1996). We rewrite the 
nonlinear Eqs. (^9|) and ( ^l|) in the form 



(38) 



where C is the linear part and TV represents the non-linear 
terms. We search for a solution of Eqs. (|3^) as a combina- 
tion of the eigenvectors obtained from the previous linear 
analysis 
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(1) 



■nil 



(2) 



(3) 



■nil 

■n± . 



■nil 



,(39) 



where tpri{t), ^vW ^'^^ 'Pvit) time dependent amplitudes 
of eigenvectors. There are always four solutions differing 
in eigenvectors multiplied by ±i or ±1. We select only so- 
lutions with physical relevance given by eigenvectors . 
We find orthonormal vectors (E^, r2^|^ , 17^^^)'^^'^'^' in order 
that 



0) 



= 1 



v± . 



■nil 



where j,k — 1,2,3. The orthonormal vectors are 



(40) 




,(3) 



AA,(E^,r!^,i,17^,)^'M, (44) 



where we denote lo^^^ = wij^max. 

The Eq. (Q) is decoupled from the others and its solu- 
tion is the decrease of the initial value of (j). Eqs. (^) and 
( ^3| ) are coupled through the linear and nonlinear terms. 
The interaction through the linear terms is weak, since 
the coupled linear terms have smaller amplitudes com- 
pared to linear terms and they decrease with time due to 
their dependence on the time derivatives of the eigenvec- 
tors, which are very small in the later stages of the shell 
evolution. The coupling through the nonlinear terms leads 
to the terms of the third and higher orders, which can be 
neglected with respect to quadratic terms. Furthermore, 
the solution of the Eq. (^2|) has a decreasing character, 
because the first term on the right side, which includes 
the "stable" U!^^\ dominates. The Eq. ( ^ ) is the most in- 
teresting one, because it has w'^' in the first linear term, 
and only the w^^^ can be imaginary negative, which has 
meaning of instability. The explicit form of the Eg. (Esf) is 



—tL- 



(E^,0^,,,r!^j(3) ^ (0^0,1) 



J7So(tj(2)-w(i)) ' (t^(2)-w(i)) 



(41) 



We insert ansatz (3£) into Eq. (pSl), multiply it by the 
orthonormal vectors (41) and obtain a set of equations for 
amplitudes i^r,{t), S.r,(t) and 




^l/^o" ^ (E o o 

I 111 I ' V^^' "'7|| ' "'^ 



,{Er„nr„,nr,^f^\^r, (42) 




(2) 



,(E^,0^ii,l]^J^'M (43) 



>(2) 



,VR-V' 



(|r,|Eo + |r,|Eo)(c^(i)-*2^) + |r7|Eo 



|j7|Eo(w(i) -cj(2)) 



;(2)) 



|r,|Eo (c.W-c.(2))'^ 



i\T] - r/'lEo {ipn^n' + S.rj-r,' 



- ^^'^ ) Crj-r,' 



i?Eo - cj(2)) ^ 



(45) 



l»7-»7'|So ('0i7-'7' +^77-77' 



i2 uj 

R 



(1) 



V 

i2 iv 

R 



(2) 



\V'\\V\ \v'^\\ri\ 



V 



E 



V 



i2 UJ"- 

R 



(2) f 

(^-^7'±,T7') 



-n 7| h (Pry-T?' 

|J7-T7'| 



V 

^2--u; 



V 

r' 
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Fig. 3. The evolution of the maximum perturbation of 
the surface density in the case when the initial values of 
linear and nonlinear terms of perturbation are in phase. 



,(2) 



Wi_\\v\ 



E 



,(2)' 



\v-v\\v\ 



v~v'±\\v\ 



.2^ 
R 



,(1) 



y 



ii^ri-r,' + £,ri-r,') \v'\^0 (V'r,' + Cry' 



In analogy to Fuchs ( 1996 ) we group together compo- 
nents with fully imaginary into wave-packets and take 
the wavenumber ?7max- Other components are grouped to 
the wave-packets of the same width and approximated by 
the average wave-numbers. The ^tj^^^ modes grow, the 
V'Tj.nax ^iid the 4>r),-n^^ modes descend. Other modes oscil- 
late with decrease. 

Demanding 
geometrical consideration 



77' = |t7 — Jy'l = ?7max we obtain from 



{v-v'±,v'±) = 



,1„2 
2 'max 



2 'max J 



(46) 



where r]'_^ and rj'^ are two wave- vectors inclined at angles 
60° to the wave-vector ij. It means that £^ri wave-packets 
non- linearly interact with others with wave- vectors rj'j.. 
Using Eqs. (|9|) and (|6|), the set of Eqs. (^ can be sim- 
plified to the form 



(47) 
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Fig. 4. The solution of the set of Eqs. (47) with the initial 
conditions corresponding to Fig. ^. 

where ^ = C,,, C+ = Cr)+, = frj-i the asterisk has 
meaning of the complex conjunction and 



A ^ -- 



2Eo(9V^2c2 -h7r2G2p2i?4) 

idRVc^^Q-/ + 
2E;^(97VT^2^2^5p4y 



and 

B ^ 



■*7) 



4c47 



where 



7 = 



i?2 



7r2G2E2 



(48) 



(49) 



(50) 



The set of Eqs. (^Tj) describe the time evolution of 
one triplet of the most interacting modes (?7, 77+, r/_). A 
is coming from the time derivative of the eigenvector in 
(^3|). Its amplitude is much less than the amplitude of w^^^ 
except around the time ^b, when w'^' is close to zero. The 
resulting effect is the time evolution of C,C+,C_, which 
for t > initially decrease and start to grow only with 
some delay (see Fig. 4). Nevertheless, the transformation 
of C,C+,C_ to S gives the continuous increase of S after 
ib (see Fig. 3). A is the term originating in the transfor- 
mation of E, il to C, and it is important around ib 
only. 

6.1. The numerical solution 

The set of Eqs. ( p7| ) can be solved numerically. We start at 
the time t\, which is the time when the instability begins 
(imaginary part of w^^max starts to be negative). First we 
select real and imaginary parts of all initial perturbation 
amplitudes C, which have the meaning of initial 

perturbations of the surface density and of the velocity. 
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Fig. 5. The evolution of the maximum perturbation of the 
surface density in the case when the initial values of linear 
and nonlinear terms of perturbation are in anti-phase. 



such that they correspond to Si/Eo = 0.05. The mag- 
nitude of these perturbations in physical values can be 
computed from the eigenvectors (|3l|). 

The solution is determined by parameters of two types: 
the first ones, as speed of sound c in the shell, are constant 
values, the second ones, as the radius of the shell R{t), 
the expansion velocity V{t) and its time derivative and 
the surface density So(t) and its time derivative, are func- 
tions of time. We can get them either from the analytical 
Sedov solution or from the numerical simulations 

of the expanding HI shells described by Ehlerova et al. 
( |1997D . In this paper we use the Sedov solution (Eqs. H-ph 
with following parameters: total energy Etot = 10^'^rg, 
density of ambient medium no = 1 cm~'^, average molecu- 
lar weight fi = 1.3, sound speed in the shell c = 1 kms~^. 

The time evolution of Sq and Si are presented in 
Figs. ^, |[ and |^ and the corresponding amplitudes of real 
and imaginary parts of ^,^_|_,^_ for the first two cases in 
Figs. I and I We can distinguish two situations: the linear 
and non-linear parts of the perturbation are in phase, so 
that they support each other, which is seen in Figs. I and 
^, or they are in anti-phase, so that the nonlinear contri- 
bution slows down the linear growth of perturbation, as it 
is visible in Figs. ^ and ^. Contribution of the non-linear 
terms depends on the shape of the forming fragments, i.e. 
on the value of the amplitude functions ^, ^+ , ^_ . Figs. ^- 1| 
show the extreme cases of that contribution. Intermediate 
cases, keeping the initial value of the perturbation in sur- 
face density at the same level, Si/Sq = 0.05, are given in 
Fig. 0. We can also see in Fig. || that the maximum con- 
tribution of nonlinear terms to the value of the perturbed 
surface density, at the time when Sq = Si, is 25% of 
the linear value. 

All the curves start at the instability time th- functions 
Ejj, Ttri^ and Sj,_ grow, although the absolute value of 
the appropriate amplitude functions ^, ^+ and ^_ descend 
during the short time after the instability begins. It is 
because S^, E^^ and are connected to the amplitude 
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Fig. 6. The solution of the set of Eqs. (47) with the initial 
conditions corresponding to Fig. 
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Fig. 7. The evolution of the maximum perturbation of 
the surface density for intermediate phase-shifts between 
initial values of the linear and nonlinear terms of pertur- 
bation. 



functions through the eigenvector (|3l|), whose E part is 
always growing with time. 

The surface density S at any point of the tangential 
plane may be computed at any expansion time after tb 
using eigenvectors (|3T| ) and Eq. (^8|) written for modes 
(?7, rj+, J7_). In Fig. gwe show the distribution of the sur- 
face density S and in Fig. |9| the velocity field of the surface 
flows V in the tangential plane for rymax in the case when 
the initial perturbations have linear terms in phase with 
the nonlinear terms corresponding to Figs. ^ and |^ at the 
time i = 55 Myr. Figs. 10 and |l^ give E and v in the 
tangential plane at the same time for the case when the 
linear and nonlinear terms of the initial perturbation are 
in anti-phase corresponding to Figs. ^ and|^. In the former 
case, the fragments are well defined and the density peaks 
are separated one from another with deep depressions in 
E. In the later case, there are high surface density chains 
with no distinct peaks and we cannot separate individual 
fragments. 
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Fig. 8. The distribution of S in the tangential plane at 
the time t = 55 Myr for tlic perturbation shown in Figs. || 
and|. 




Fig. 9. The velocity vectors of the surface flows v corre- 
sponding to Fig. 0. 



We measure the total mass concentrated in one of the 
well defined fragments shown in Fig. ||. The mass of this 
fragment is given in Fig. ^ as a function of time. It de- 
creases because the decrease of its size, which is propor- 



tional to An 



and increases because the accu- 



mulation of the amfent medium and surface flows. After 
th the resulting mass of the fragment decreases, since the 
influence of the size shrinking dominates. This happens 
when the magnitude of ^^^ff^ is larger then the amplitude 
of the surface flows Umax (see Fig 
of ^^^ffisx decreases with time and i 



dt 



13). The magnitude 
lax increases, and at 
t ~ 53 Myr they are equal. Since then the inflow domi- 
nates and the fragment mass growths. 

6.2. The mass spectrum of fragments 

After tb, when the first mode begins to be gravitationally 
unstable, more and more modes are unstable and the in- 
terval of instability growths. In Fig. [ij we give the values 



Fig. 10. The same as in Fig. ^ at the same time for the 
perturbation shown in Figs. |^ and ^. 




Fig. 11. The velocity vectors of the surface flows v corre- 
sponding to Fig. |lO|. 

of the fragmentation integral Ii,T]it) as defined in ( ^ ) as 
a function of time. This shows at any time the level of 
development of a fragment with given r/. 

Mass mfrag of a fragment, which is related to A = 
may be defined as 



TOfra 



A\% _ n^R'^o 



(51) 



The formation frequency of fragments corresponding to 
the wavenumber 77 is proportional to li^ri, and the surface 
of the spherical shell of radius R is able to accommodate 
R? j}? fragments of wavelength A. Therefore, the number 
diV of fragments of mass mfrag formed out of the spheri- 
cal shell of radius R in gravitational unstable modes with 
wavenumber r] from the interval (ry, 77 + Ar\) is 



diV 



(52) 



With this equation, we may derive the mass spectrum of 
fragments dN/drrifyag at some time. It is shown in Fig. 



R. Wiinsch and J. Palous; Gravitational instability of expanding shells 



9 



20000 





Fig. 12. The time evolution of total mass in a well defined pjg_ ^hc time evolution of the fragmentation integral 
fragment. j^^^^ gi^^n by Eq. (||). 
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Fig. 13. A and Umax as functions of time. 



at the time if, when Sg = ^i- The masses are between 
a few times 10^ Mq and 1.5 x 10'^ Mq with the highest 
probability peak at Wfrag — 10^ Mq. At that time is the 
radius of the shell almost Ikpc with the expansion velocity 
~ 5 kms^^, and the total mass collected in the shell is 
1.14 X 10^ Mq. 

The decreasing part of the mass spectrum can be ap- 
proximated as a power law dA^/drntiag ~ '^^frag- 
of the this part of dA^/drrifrag gives a = —1.4, which is 
close to the ob served mass spectrum of GMC in the Milky 
Way: Combes( |l99l| ) gives a = -1.5. NANTEN survey of 
the CO emission of the LMC (Fukui, 2001) gives steeper 
slope of a = —1.9, which may be explained in the con- 
nection to higher level of random velocities in the LMC 
compared to the Milky Way resulting in the deficit of high 
mass clouds. 



6.3. The time of fragmentation 

The evolution of the maximum perturbation of the sur- 
face density can be used to determine the fragmentation 
time ti of the shell. Because at advanced stages the value 




Fig. 15. The mass spectrum of fragments. The straight 
line is the power law fit of the decreasing part of the spec- 
trum m~ 
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Fig. 16. Dependence of the fragmentation time on the ini- 
tial perturbation of the surface density. The vertical line 
gives the time when the value of /f,j;„,^x = 1. 



10 



R. Wiinsch and J. Palous; Gravitational instability of expanding shells 



of the maximum perturbation rises steeply with the time 
(see e.g. Fig. ||), we define tf as the time, when maximum 
perturbation of the surface density is equal to the unper- 
turbed value: Ei(if) = So(tf). Using the fragmentation 
integral If^ri we may compare the development level of dif- 
ferent fragments at tf (see Fig. ^4|). We can say that the 
most frequent fragments are also the most developed, the 
more massive form only later. 

tf depends on the initial conditions of the set of Eqs. 
(^). They correspond to the initial perturbation of the 
surface density. We can set them to the value typical for 
the inhomogeneities in the clumpy interstellar medium 
(10^9 - 10^0 cm-2), which is at t^: 0.01 - 0.2 x Eq- 

The dependence of tf on the value of the initial pertur- 
bation, Si(tb), is shown in Fig. Fragments form since 
45 Myr, for the largest perturbations, to 70 Myr, for the 
smallest perturbations. The spread in tf for given £i(tb) 
is connected to the different shape of the perturbation 
as shown in Fig. ^. This time may be compared to the 
fragmentation time tf^^i obtained for rymax from the linear 
analysis defined as a time when If^^i — 1- 

7. Conclusions 

We evaluate the time evolution of perturbations on the 
surface of an expanding shell. We complement the linear 
analysis of the gravitational fragmentation process with 
the inclusion of nonlinear terms, and we compute the time 
evolution of fragments after the time when the shell starts 
to be unstable. Some initial perturbations develop into 
well separated fragments and we estimate the time evo- 
lution of the mass of a fragment, the mass spectrum of 
fragments, and the spread in their formation time. The 
computed mass spectrum is close to the observed mass 
distribution of GMC in the Milky Way, but slightly flat- 
ter than the mass spectrum of molecular clouds observed 
in the LMC. This may be related to higher level of random 
motions in the LMC compared to the Milky Way, which 
restricts the formation of late time massive fragments and 
steepen the resulting mass spectrum. Also interesting is 
that the more massive fragments form at later times of the 
shell evolution than the less massive fragments. The for- 
mation time depends on the value of the initial perturba- 
tion: tf = 45 — 70 Myr. Large density fluctuations shorten 
this time and thus in the disturbed ISM with large density 
fluctuations the fragments form sooner than in quiet and 
smooth ISM where the density fluctuations are small. 
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